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We present new results on chiral symmetry breaking in nearly conformal gauge theories with 
fermions in the fundamental representation of the SU(3) color gauge group. The number of 
fermion flavors is varied in an extended range below the conformal window with chiral sym- 
metry breaking (^SB) for all flavors between Nf = 4 and Nf = 12. To identify jSB we apply 
Q\ several methods which include, within the framework of chiral perturbation theory, the analysis 

of the Goldstone spectrum in the p-regime and the spectrum of the fermion Dirac operator with 
eigenvalue distributions of random matrix theory in the e-regime. Chiral condensate enhance- 
ment is observed with increasing Nf when the electroweak symmetry breaking scale F is held 
fixed in technicolor language. Important finite-volume consistency checks from the theoretical 
understanding of the SU(Nf) rotator spectrum of the 5 -regime are discussed. We also consider 
these gauge theories at Nf — 16 inside the conformal window. The importance of understanding 
finite volume, zero momentum gauge field dynamics inside the conformal window is pointed out. 
Staggered lattice fermions with supressed taste breaking are used throughout the simulations. 
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1. Introduction 

Our goal in this work to identify chiral symmetry breaking (#SB) below the conformal win- 
dow of strongly interacting gauge theories requires the application and testing of several methods 
in finite volumes. These include the analysis of the Goldstone spectrum in the p-regime and the 
spectrum of the fermion Dirac operator with eigenvalue distributions of Random Matrix Theory 
(RMT) in the £-regime, within the framework of chiral perturbation theory (#PT). Some critical 
consistency checks from the theoretical understanding of the SU(Nf) rotator spectrum of the 8- 
regime will be also discussed. We report new results with Nf = 4, 8,9, 12 flavors with #SB below 
the conformal window for fermions in the fundamental representation of the SU(3) color gauge 
group. As Nf is increased, chiral condensate enhancement is observed when the electroweak sym- 
metry breaking scale F is held fixed in technicolor language. We also discuss the theory inside the 
conformal window. The importance of understanding finite volume, zero momentum gauge field 
dynamics inside the conformal window is pointed out and illustrated at Nf =16. Much of this work 
is an extension of our pre-conference publication [1] where we did not report our Nf = 12 results. 
In a forthcoming publication [2] more details will be provided on the analysis and results presented 
here. Our work on the running coupling is presented separately [3] . 

It is an intriguing possibility that new physics beyond the Standard Model might take the 
form of some new strongly-interacting gauge theory building on the original technicolor idea [4, 
5, 6]. This approach has lately been revived by new explorations of the multi-dimensional theory 
space of nearly conformal gauge theories [7, 8, 9, 10]. Model building of a strongly interacting 
electroweak sector requires the knowledge of the phase diagram of nearly conformal gauge theories 
as the number of colors N c , number of fermion flavors Nf, and the fermion representation R of the 
technicolor group are varied in theory space. For fixed iV c and R the theory is in the chirally broken 
phase for low Nf, and asymptotic freedom is maintained with a negative /3 function. On the other 
hand, if Nf is large enough, the /3 function is positive for all couplings, and the theory is trivial. 
There is some range of Nf for which the j8 function might have a non-trivial zero, an infrared 
fixed point, where the theory is in fact conformal [13, 14]. This method has been refined by 
estimating the critical value of Nf, above which spontaneous chiral symmetry breaking no longer 
occurs [15, 16, 17]. 

Interesting models require the theory to be very close to, but below, the conformal window, 
with a running coupling which is almost constant over a large energy range [18, 19, 20, 21, 22, 23]. 
The nonperturbative knowledge of the critical Nj r " separating the two phases is essential and this 
has generated much interest and many new lattice studies [1, 24, 25, 26, 27, 28, 29, 30, 31, 33, 
32, 34, 35, 36, 37, 38, 39, 40, 41, 42, 43, 44, 45, 46, 47, 48, 49, 50, 51, 52, 53, 54, 55, 56]. To 
provide theoretical framework for the analysis of simulation results, we review first a series of tests 
expected to hold in the setting of #PT in finite volume and in the infinite volume limit. 

2. Chiral symmetry breaking below the conformal window 

We will identify in lattice simulations the chirally broken phases with Nf ■ = 4, 8,9, 12 flavors 
of staggered fermions in the fundamental SU(3) color representation using finite volume analysis. 
We deploy staggered fermions with exponential (stout) smearing [57] in the lattice action to reduce 
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well-known cutoff effects with taste breaking in the Goldstone spectrum [58]. The presence of taste 
breaking requires a brief explanation of how we apply staggered #PT in our analysis. The important 
work of Lee, Sharpe, Aubin and Bernard [59, 60, 61] is closely followed in the discussion. 

2.1 Staggered chiral perturbation theory 

Starting with N f = 4 [59], the spontaneous breakdown of SU{A) L x SU{A) R to vector SU{4) 
gives rise to 15 Goldstone and pseudo-Goldstone modes, described by fields fa that can be orga- 
nized into an SU (4) matrix 

E(x)=e X p(V^-), <j> = J £fa,T a . (2.1) 



a=\ 



F is the Goldstone decay constant in the chiral limit and the normalization T a = { ^ , i£ps , , ^5 } 
is used for the flavor generators. The leading-order chiral Lagrangian is given by 

J^ 4) = ^-Jr(d^d^) - l -Bm q F 2 J r (L + L Jf ) , (2.2) 

with the fundamental parameters F and B measured on the technicolor scale Ajc which replaced 
Aqcd hi the new theory. Expanding the chiral Lagrangian in powers of <p one finds 15 degenerate 
Goldstone pions with masses given by 

M 2 n = 2Bm q [1 + 0K/A T c)] . (2.3) 

The leading-order term is the tree-level result while the corrections come from loop diagrams and 
from higher-order terms in the chiral Lagrangian. The addition of a 2 «Sf^ 6 ' breaks chiral symmetry 
and lifts the degeneracy of the Goldstone pions. Adding correction terms to Eq. (2.3) yields 

M 2 n = C{T a ) ■ a 2 h\ c + 2Bm q [l + 0(m q /A TC ) + 0{a 2 A 2 c )} (2.4) 

where the representation dependent C{T a ) is a constant of order unity. Contributions proportional 
to a 2 are due to and lead to massive pseudo-Goldstone pions even in the m q — > chiral limit, 
except for the Goldstone pion with flavor ^5 which remains massless because the U(1)a symmetry 
is protected. 

Lee and Sharpe observe that the part of without derivatives, defining the potential , 
is invariant under flavor SO(4) transformations and gives rise to the a 2 term in M 2 .. Terms in 
Jfj,^ involving derivatives break SO(4) further down to the lattice symmetry group and give rise 
to non-leading terms proportional to a 2 m and a 4 . The taste breaking potential is given by 



r z (6) = CiTrfeX&Zt) + C 2 \ [Tr(I 2 ) - Tr(feZ^E) +h 



+ C 3 \ £ [Tr(^ v I) + h.c] + C 4 ^ £ [Tr(£ v5 £kv£) + h.c] 

V V 

+ C 5 ^ [Tr^vtf) - Tr^vsE&vtf)] + Q, £ Tr^I^) . (2.5) 



v n<v 
The six unknown coefficients Q are all of size A^ c . 



3 



Chiral symmetry breaking in nearly conformal gauge theories 



In the continuum, the Goldstone pions form a 15-plet of flavor SU (4) and are degenerate. On 
the lattice, states are classified by the symmetries of the transfer matrix, and the pseudo-Goldstone 
pions fall into 7 irreducible representations: four 3 -dimensional representations with flavors £,-5, 
and ^4, and three 1 -dimensional representations with flavors £4, ^45 and £5. 

Close to both the chiral and continuum limits, the pseudo-Goldstone masses are given by 

M K (T a ) 2 = 2Bm q + a 2 A{T a ) + 0{a 2 m q ) + <9(a 4 ) , (2.6) 

with A(T a ) ~ A^ c arising from "fjjP . Since respects flavor SO(4), the 15 Goldstone particles 
fall into SO(4) representations: 

A(&) = 0, A(^) = A ( Cl + C 2 + C 3 + 3C 4 + C 5 + 3C 6 ) , 

A(^ 5 ) = ^(Ci+C 2 + 3C3+C 4 -C 5 +3C 6 ), A(^ V ) = A(2C 3 + 2C 4 +4C 6 ). (2.7) 

In the chiral limit at finite lattice spacing, the lattice irreducible representations with flavors and 
(§4 are degenerate, those with flavors ^5 and ^45, and those with flavors and <^ 4 are degenerate 
as well. No predictions can be made for the ordering, splittings, or even the signs of the mass 
shifts. Our simulations indicate that they are all positive with the exponentially smeared staggered 
action we use, making the existence of an Aoki phase [59] unlikely. The method of [59] has 
been generalized to the Nf > 4 case [60, 61] which we adopted in our calculations with help from 
Bernard and Sharpe. The procedure cannot be reviewed here but it will be used in the interpretation 
of our Nf = 8 simulations. 

2.2 Finite volume analysis in the p-regime 

Three different regimes can be selected in simulations to identify the chirally broken phase 
from finite volume spectra and correlators. For a lattice size L? s x L t in euclidean space and in the 
limit L t S> L s , the conditions F n L s > 1 and M n L s > 1 select the the p-regime, in analogy with low 
momentum counting [62, 63] . 

For arbitrary Nf, in the continuum and in infinite volume, the one-loop chiral collections to 
M K and F n of the degenerate Goldstone pions are given by 



M\=M 1 



(2.8) 



ln( A3 
~ 8n 2 N f F 2 \~M 

, N f M 2 , /A 4 \ 



(2.9) 



where M 2 = 2B • m q and F, B, A3 , A4 are four fundamental parameters of the chiral Lagrangian, and 
the small quark mass m q explicitly breaks the symmetry [64]. The chiral parameters F,B appear 
in the leading part of the Lagrangian in Eq. (2.2), while A3,A4 enter in next order. There is the 
well-known GMOR relation ~L com i = BF 2 in the m q — > limit for the chiral condensate per unit 
flavor [65]. It is important to note that the one-loop correction to the pion coupling constant F n is 
enhanced by a factor Nf compared to M 2 . The chiral expansion for large Nf will break down for F K 
much faster for a given M n jF n ratio. The NNLO terms have been recently calculated [66] showing 
potentially dangerous N 2 corrections to Eqs. (2.8,2.9). 



4 



Chiral symmetry breaking in nearly conformal gauge theories 



The finite volume corrections to M n and F n are given in the p-regime by 



M n {L s ,r]) =M n 



F n (L s ,r}) = F n 



1 + 



1 



1 



M 



2N f \6k 2 F 2 



N f M 2 



(2.10) 



(2.11) 



2 16k 2 F 2 

where gi {X , r\) describes the finite volume connections with A = M ■ L s and aspect ratio v\ = L t /L s . 
The form of g\{X, f]) is a complicated infinite sum which contains Bessel functions and requires 
numerical evaluation [63]. Eqs. (2.8-2.1 1) provide the foundation of the p-regime fits in our simu- 
lations. 



chiral p-regime 




Figure 1: Schematic plot of the regions in which the three low energy chiral expansions are valid. The 
vertical axis shows the finite temperature scale (euclidean time in the path integral) which probes the rotator 
dynamics of the 5-regime and the e-regime. The first two low lying rotator levels are also shown on the 
vertical axis for the simple case of Nf = 2. The fourfold degenerate lowest rotator excitation at m q = 
will split into an isotriplet state (lowest energy level), which evolves into the p-regime pion as m q increases, 
and into an isosinglet state representing a multi-pion state in the p-regime. Higher rotator excitations have 
similar interpretations. 

2.3 5-regime and £-regime 

At fixed L s and in cylindrical geometry L,/L s S> 1, a crossover occurs from the p-regime to 
the 5 -regime when m q — > 0, as shown in Fig. 1. The dynamics is dominated by the rotator states 
of the chiral condensate in this limit [67] which is characterized by the conditions FL S > 1 and 
ML S <C 1. The densely spaced rotator spectrum scales with gaps of the order ~ \/F 2 L? s , and at 
m q = the chiral symmetry is apparently restored. However, the rotator spectrum, even at m q = 
in the finite volume, will signal that the infinite system is in the chirally broken phase for the 
particular parameter set of the Lagrangian. This is often misunderstood in the interpretation of 
lattice simulations. Measuring finite energy levels with pion quantum numbers at fixed L s in the 
m q — ► limit is not a signal for chiral symmetry restoration of the infinite system [41]. 
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0.16 




Figure 2: The crossover from the p-regime to the 5 -regime is shown for the n and states at Nf = 4. 

If L t ~ L s under the conditions FL S > 1 and ML S <C 1 , the system will be driven into the 
£-regime which can be viewed as the high temperature limit of the 5-regime quantum rotator. 
Although the 5-regime and £-regime have an overlapping region, there is an important difference 
in their dynamics. In the 5-regime of the quantum rotator, the mode of the pion field U (x) with zero 
spatial momentum dominates with time-dependent quantum dynamics. The £ -regime is dominated 
by the four-dimensional zero momentum mode of the chiral Lagrangian. 

We report simulation results of all three regimes in the chirally broken phase of the technicolor 
models we investigate. The analysis of the three regimes complement each other and provide 
cross-checks for the correct identification of the phases. First, we will probe Eqs. (2.8-2.11) in 
the p-regime, and follow with the study of Dirac spectra and RMT eigenvalue distributions in the 
£-regime. The spectrum in the c> -regime is used as a signal to monitor p-regime spectra as m q 
decreases. Fig. 2 is an illustrative example of this crossover in our simulations. It is important 
to note that the energy levels in the chiral limit do not match the rotator spectrum at the small 
F ■ L s values of the simulations. This squeezing with F • L s not large enough for undistorted, finite 
volume, chiral behavior in the p-regime, £-regime, and 5-regime, will be further discussed in our 
p-regime simulations presented next. We will also describe some methods to put this squeezing 
into a more quantitative context. 

3. Goldstone spectra and #SB from simulations at Nf = 4 in the p-regime 

In this section we describe in some detail the methods we use for successfully testing chiral 
symmetry breaking. Our tests in the p-regime have two major components. The primary test is to 
identify the pseudo-Goldstone spectrum of the staggered formulation with evidence for recovery 
from taste symmetry breaking close to the continuum limit. The secondary test is to probe chiral 
loop corrections to the tree-level behavior of M\ and F n as the fermion mass a ■ m q is varied at 
fixed gauge coupling /3. The evidence we find for chiral symmetry breaking at Nf = 4, 8,9, 12 is 
common to all flavors. Limitations and ambiguities identified at Nf = 4 for future improvements are 
expected to be more pronounced with increasing Nf. Results for each flavor we have simulated in 
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Figure 3: The Goldstone spectrum and chiral fits are shown for simulations with lattice size 24 3 x 32. The 
left column shows the pseudo-Goldstone spectrum with decreasing taste breaking as the gauge coupling is 
varied from /3 = 3.5 to /3 = 3.7. The middle value at /3 =3.6 was chosen for chiral fits which are shown 
in the right column. The top right figure with fitting range a ■ m q — 0.008 — 0.025 shows the NLO chiral 
fit to M\jm q which approaches IB in the chiral limit. Data points below m q = 0.008 are not in the chiral 
p-regime and not used in the fitting procedure. The middle figure on the right is the NLO chiral fit to F % in 
the range a ■ m q = 0.008 — 0.02. The bottom right figure is the linear fit to the chiral condensate with fitting 
range a-m q = 0.015 — 0.025. The physical fit parameters B,F, A3, A4 are discussed in the text. 

the p-regime are presented in separate sections beginning here with general discussion and Nf = 4 
results. 

We have used the tree-level Symanzik-improved gauge action for all simulations. The con- 
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ventional /3 = 6/g 2 lattice gauge coupling is defined as the overall factor in front of the three 
well-known terms of the lattice action. The link variables in the staggered fermion matrix were 
exponentially smeared with six stout steps at Nf = 4 and the RHMC algorithm was deployed in all 
runs. The results shown in Fig. 3 are from the p-regime of #SB with the conditions M n ■ L s » 1 and 
F n ■ L s ~ 1 when the chiral condensate begins to follow the expected behavior of infinite- volume 
chiral perturbation theory from Eqs. (2.8,2.9) in next-to-leading order (NLO) with calculable fi- 
nite volume corrections from Eqs. (2.10,2.1 1) which are negligible at L s = 24. We have empirical 
evidence that the M n and F % data points are free of finite volume corrections in practically the en- 
tire fitting range of the fermion masses we use at L s = 24 so that the negligible corrections from 
Eqs. (2. 10,2. 1 1) can be ignored. 

Within some finite volume limitations, which we will address, the Nf = 4 simulations work 
in the p-regime as expected. The left column of Fig. 3 shows that the pseudo-Goldstone spectrum 
clearly remains separated from the hadronic scale of the p -meson as j3 is varied. Moving towards 
the continuum limit with increasing jS , we see the split pseudo-Goldstone spectrum collapsing into 
the degenerate continuum pion spectrum. The true Goldstone pion whose mass will vanish in the 
a-m q = limit at fixed lattice spacing and two additional split states with small residual masses at a ■ 
m q = are shown to illustrate the trend, a 4 A is the measure of the small taste breaking in quadratic 
mass splitting as measured in lattice units. The origin of the splittings and the quantum numbers 
were discussed in Section 2 as shown in Eq. (2.7). The spectrum is approximately parallel as the 
bare fermion mass a • m q is varied at fixed lattice spacing and the gaps appear to be equally spaced 
to a good approximation, consistent with earlier observations in QCD where the C4 term seems 
to dominate staggered taste breaking for two light flavors with equally spaced pseudo-Goldstone 
levels [59]. We selected /3 = 3.6 for testing #PT of finite volume Goldstone dynamics in the p- 
regime. This choice with small taste breaking is close to the continuum limit without excessive 
squeeze on the important product F ■ L s which in an ideal simulation of the p-regime should be 
large (F is the chiral limit of F n as a ■ m q — ► at fixed lattice spacing). 

The simultaneous chiral fit of M 2 /m q and F K based on Eqs. (2.8-2.11) is shown in Fig. 3 
where chiral loops correct the tree-level values of M 2 /m q = IB and F n . In the fitting range a-m q = 
0.008 - 0.025 applied to M\jm q we observe small corrections to the tree-level value of IB which 
keeps the fit well within the range of one-loop jYY. In the fitting range a-m q = 0.008 — 0.02 the 
F n data are about a factor of two larger than F which indicates how the one-loop fit is being pushed 
to its limits. Without loop correction F K would not change from its fitted value of a ■ F = 0.033(4) 
in the chiral limit at fixed lattice spacing. The fitted value of B is a B = 1.76(7) in lattice units 
and Mp/F = 13(1) in the chiral limit (the linear fit of M p = c + d • m q is used at all Nf values to 
determine M p (m q = 0)). The fitted value of B/F = 53(6) indicates significant enhancement of the 
chiral condensate from its Nf = 2 value [53, 68]. In our simultaneous fits we get A3 = 0.37(5) and 
A4 = 0.51(1) which set the chiral couplings in the NLO chiral Lagrangian. 

The chiral condensate (yy) summed over all flavors is dominated by the linear term in m q 
from UV contributions. The linear fit gives (yy) = 0.0191 (4) in the chiral limit which differs from 
the GMOR relation of (yy) = 4F 2 B by about a factor of two with 4F 2 B = 0.008(2) fitted. There 
are several sources of this disagreement. The chiral log in (yy) will bring further down the true 
fitted value in the chiral limit. Our volumes are not large enough yet to attempt a sensible chiral log 
fit to the condensate at small a ■ m q values. Finite volume squeezing effects distort the consistency of 
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the results in our limited range of simulation volumes. The choice of fitting method to Eqs. (2.8,2.9) 
can also have some effect on the results. On the right-hand sides of the equations, the variable pair 
(M, F) in the chiral logs can be replaced with the pair (M n , F n ) which is equivalent to a partial 
resummation [68]. This will be reported in our more detailed forthcoming journal publication [2]. 

Finite volume limitations when measured in F units have the biggest effect on our chiral anal- 
ysis. The value F ■ L s « 0.8 is not expected to protect against significant finite-volume squeezing 
effects for even the largest spatial size L s = 24 used in the chiral fits at Nf = 4. Larger than opti- 
mal NLO corrections in the chiral fits of F K and finite- volume squeezing effects are closely related 
concerns. Simulations on larger lattices would increase F ■ L s and allow us to drop back in a ■ m q 
into a more comfortable range with smaller NLO chiral corrections for F n . 

Finite volume corrections to the rotator spectrum in the 5 -regime set some quantitative mea- 
sure of squeezing effects on the chiral analysis. The connection is made by observing that the 
pion spectrum in the p-regime can be viewed at fixed spatial volume h\ as the adiabatic evolu- 
tion from the energy levels of the rotator spectrum of the 5 -regime as illustrated schematically 
in Fig. 2 for the lowest Nf = 2 rotator levels. The rotator spectrum for Nf = 2 is given by 
Ei = 4aZ(/ + 2), with / = 0,1,2,..., where the moment of inertia is calculated in NLO [69, 70] 
as = F 2 L](\ + C( ^g 2) + 0(^4 )). The value of C(N f = 2) is known to be 0.45 and is expected 

to grow with Nf. At F • L s m 0.8 the correction is 70% for Nf = 2 and probably considerably larger 

Ni-l 

for Nf = 4. The leading-order rotator gap for arbitrary Nf is given by E\—Eq= 1 2 3 but the 

l\fr L s 

coefficient C(Nf) is an important missing piece in the analysis. Were we to continue the p-regime 
Goldstone spectrum at Nf = 4, L s = 24, and j8 = 3.6 to the 5 -regime adiabatically, the small value 
of F ■ L s would not allow us to get a reliable estimate of F based on the chiral rotator spectrum 
with the collapse of the adiabatic approximation. This is a quantitative warning sign of the need for 
considerably larger spatial volume for robust p-regime results to determine F in the chiral fitting 
procedure. In fact, we are going to fit M\jm q in the chiral analysis for Nf = 8,9, 12 with better 
controlled NLO loop corrections, but F K will not be fitted. For larger Nf, a reliable simultaneous 
fit requires substantially larger volumes than are realistic with our current resources. 

In summary, the Nf = 4 system passed both tests in the chirally broken phase and shows signif- 
icant enhancement of the chiral condensate when measured in units of the electroweak symmetry 
breaking scale set by F. This is a relevant effect to monitor for fermion mass generation in extended 
technicolor applications as we begin to approach the conformal window [53]. 

4. Goldstone spectrum and #SB from simulations at Nf = 8 in the p-regime 

As we move to the Nf = 8 p-regime simulations, we can clearly identify the p-regime of the 
chirally broken phase as summarized in Fig. 4. The same lattice action and algorithm was used 
for the Nf = 8 p-regime simulations as introduced earlier for Nf = 4. We can clearly identify the 
pseudo-Goldstone spectrum which is separated from the technicolor scale of the p -meson. Moving 
towards the continuum limit we observe at /3 = 1 .4 the split pion spectrum collapsing toward the 
true Goldstone pion with a new distinct feature. The true Goldstone pion, whose mass will vanish 
in the a • m q = limit at fixed lattice spacing, and two additional split pseudo-Goldstone states 
appear with considerably different slopes in Fig. 4 as m q increases. For small a ■ m q we find the 



9 



Chiral symmetry breaking in nearly conformal gauge theories 




Figure 4: The Goldstone spectrum and chiral fits are shown for Nf = B simulations with lattice size 24 3 x 32. 
The left column shows the pseudo-Goldstone spectrum with decreasing taste breaking as the gauge coupling 
is varied from /3 = 1 .0 to j5 — 1 .8. The middle value at j3 — 1 .4 was chosen in the top right figure with fitting 
range a ■ m q = 0.015 — 0.03 of the NLO chiral fit to M^/m q which approaches 2B in the chiral limit. The 
middle figure on the right shows the F K data with no NLO fit far away from the chiral limit. The bottom 
right figure is the linear fit to the chiral condensate with fitting range a-m q = 0.02 — 0.04. The physical fit 
parameters B,F,Ai are discussed in the text. 



pseudo-Goldstone spectrum collapsed at fixed gauge coupling. Apparently the NLO operators, 
the last two terms in Eq. (2.6), have a stronger effect on the spectra relative to leading-order taste 
breaking operators, the generalization of those from Nf = 4 to Nf = 8 as discussed in Section 2. 
This somewhat unexpected and unexplained trend is observed for Nf > 8 as well. 
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We analyzed the #SB pattern within staggered perturbation theory in its generalized form 
beyond four flavors [60, 61]. The simultaneous chiral fit of M 2 /m q and F n based on Eqs. (2.8-2.1 1) 
cannot be done at Nf = 8 within the reach of the largest lattice sizes we deploy since the value of 
F • L s is too small even at L s = 24, for gauge couplings where taste breaking drops to an acceptable 
level. The chiral fit of B,F,A^ to M 2 /m q , shown at the top right of Fig. 4, is based on Eq. (2.8) only 
since the F K data points are outside the convergence range of the chiral expansion. Much larger 
lattices are required to drop down in m q to the region where the simultaneous fit could be made, 
while maintaining some control over finite volume and taste breaking effects. The finite volume 
corrections were negligible in the fitted a • m q range and Eqs. (2. 10,2. 1 1) were not needed. 

At j8 = 1 .4 the fitted value of B is a ■ B = 2.6(3) in lattice units with a-F = 0.0166(9) and 
a • A 3 = 0.48(5) also fitted. The fitted p-mass in the chiral limit is a -M p = 0.27(2) with M p /F = 
17(1). The fitted value of B/F = 158(17) is not very reliable but indicates significant enhancement 
of the chiral condensate from its Nf = 4 value without including renormalization scale effects. 
The chiral condensate (yy) summed over all flavors is dominated by the linear term in m q from 
UV contributions. The linear fit gives (yy) = 0.012(1) in the chiral limit which differs from the 
GMOR relation of (yy) = SF 2 B by about a factor of two with SF 2 B = 0.0058(8) fitted. There are 
several sources of this disagreement which were addressed for the Nf = 4 case earlier. The chiral 
log in (yy) will bring further down the true fitted value in the chiral limit. Our volumes are not 
large enough yet to attempt a sensible chiral log fit to the condensate at small a ■ m q values. Finite 
volume squeezing effects distort the consistency of the results in our limited range of simulation 
volumes. Similar observations should also be noted when the RMT analysis is applied in the e- 
regime. 

5. Goldstone spectrum and #SB from simulations at Nf = 9 in the p-regime 

We had two motivations for the Nf = 9 simulation reported here. We wanted to see whether 
the rooting procedure (being applied in our project with two fermions in the sextet representation) 
will present some unexpected changes in the analysis and we were also looking for the continued 
trends in the #SB pattern. We could not find any noticeable effect from the rooting procedure and 
the symmetry breaking pattern was consistent with the Nf = 8 simulations. 

As shown in Fig. 5 the Goldstone spectrum is still clearly separated from the technicolor scale 
of the p -meson. The true Goldstone pion and two additional split pseudo-Goldstone states are 
shown again with different slopes as a • m q increases. The trends and the underlying explanation 
are similar to the Nf = 8 case. The chiral fit to M 2 /m q is shown based on Eq. (2.8) only since 
the F % data points are outside the convergence range of the chiral expansion. At j3 = 2.0 the fitted 
value of B is a ■ B = 2.8(4) in lattice units with a-F = 0.017(2) and a • A 3 = 0.48(9) also fitted. 
The fitted p-mass in the chiral limit is a -M p = 0.233(3) with M p /F = 14(1). The fitted value of 
B/F = 166(32) is not very reliable but comparable to the enhancement of the chiral condensate 
found at Nf = 8 without including renormalization scale effects. Again, at fixed lattice spacing, the 
small chiral condensate (yy) summed over all flavors is dominated by the linear term in m q from 
UV contributions. The linear fit gives (yy) = 0.0045(7) in the chiral limit which differs from 
the GMOR relation of (yy) = 9F 2 B by about a factor of two with 9F 2 B = 0.007(2) fitted. Open 
issues in the systematics are similar to the Nf = 8 case. 
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Figure 5: The pseudo-Goldstone spectrum and chiral fits are shown for Nf = 9 simulations with lattice size 
24 3 x 32. The left column shows the pseudo-Goldstone spectrum with decreasing taste breaking as the gauge 
coupling is varied from j3 = 1.6 to j3 = 2.4. Although the bottom figure on the left at j3 = 2.4 illustrates the 
continued restoration of taste symmetry, the volume is too small for the Goldstone spectrum. The middle 
value at j3 = 2.0 was chosen in the top right figure with fitting range a ■ m q = 0.015 — 0.03 of the NLO chiral 
fit to M\jm q which approaches 2B in the chiral limit. The middle figure on the right shows the F n data with 
no NLO fit far away from the chiral limit. The bottom right figure is the linear fit to the chiral condensate 
with fitting range a ■ m q = 0.02 — 0.04. The physical fit parameters B,F,Aj, are discussed in the text. Four 
stout steps were used in all Nf — 9 simulations. 

6. Goldstone spectrum and ^SB from simulations at Nf = 12 in the p-regime 

Finally we move to the controversial Nf = 12 case. We find here a similar chiral symmetry 
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Figure 6: The pseudo-Goldstone spectrum and chiral fits are shown for Nf =12 simulations with lattice size 
24 3 x 32 and 32 4 . The left column shows the pseudo-Goldstone spectrum with decreasing taste breaking as 
the gauge coupling is varied from j3 = 2.0 to j3 = 2.4. Although the bottom figure on the left at j3 = 2.4 
illustrates the continued restoration of taste symmetry, the volume is too small for the Goldstone spectrum. 
The middle value at j3 = 2.2 was chosen in the top right figure with fitting range a ■ m q = 0.015 — 0.035 of 
the NLO chiral fit to M\jm q which approaches 2B in the chiral limit. The middle figure on the right shows 
the F K data with no NLO fit far away from the chiral limit. The bottom right figure, with its additional 
features discussed in the text, is the linear fit to the chiral condensate with fitting range a ■ m q = 0.02 — 0.04. 
The physical fit parameters B,F,A^ are discussed in the text. Two stout steps were used in all Nf = 12 
simulations. 
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breaking pattern as we found in the Nf = 8,9 cases with increased concerns about all the caveats 
presented before. The Goldstone spectrum remains separated from the technicolor scale of the 
p -meson. The true Goldstone pion and two additional split pseudo-Goldstone states are shown 
again in Fig. 6 with different slopes as a ■ m q increases. The trends and the underlying explanation 
are similar to the Nf = S,9 cases. The chiral fit to M 2 /m q shown at the top right side of Fig. 6 
is based on Eq. (2.8) only since the F K data points are outside the convergence range of the chiral 
expansion. At j8 = 2.2 the fitted value of B is a B = 2.7(2) in lattice units with a-F = 0.0120(1) 
and a ■ A3 = 0.50(3) also fitted. The fitted p-mass in the chiral limit is a -M p = 0.115(15) from 
a-m q = 0.025-0.045 withM p /F = 10(1). The fitted value of B/F = 223(17) is not very reliable 
but consistent with the enhancement of the chiral condensate found at Nf ■ = 8,9 without including 
renormalization scale effects. Again, at fixed lattice spacing, the small chiral condensate (yy) 
summed over all flavors is dominated by the linear term in m q from UV contributions. The linear 
fit gives (yy) = 0.0033(13) in the chiral limit which came out unexpectedly close the GMOR 
relation of (yy) = 12F 2 B with \2F 2 B = 0.0046(4) fitted. Issues and concerns in the systematics 
are similar to the Nf = 8,9 cases. 

In summary, we have shown strong evidence that according to p-regime tests the Nf = 4, 8, 9, 12 
systems all exhibit in the Goldstone and hadron spectra broken chiral symmetry close to the contin- 
uum limit. There are some important features of the Nf = 12 analysis which suggest that the model 
is not only in the #SB phase but also close to slow walking of the renormalized gauge coupling. 
The bottom right of Fig. 6 shows the crossover in the chiral condensate from strong coupling to 
the weak coupling regime in the relevant range of m q . In combination with the nearly degenerate 
Goldstone spectrum we find it quite suggestive that around /3 = 2.2 we are close to continuum be- 
havior. In addition, we observe that the fitted value of the p-mass in the chiral limit hardly changes 
in this region as the gauge coupling is varied (at /3 = 2.0 we fit a -M p = 0.123(10)). If confirmed on 
larger lattices, this could be a first hint of a slowly changing gauge coupling close to the conformal 
window. Currently we are investigating the important Nf = 12 model on larger lattices to probe the 
possible influence of unwanted squeezing effects on the spectra. This should also clarify the mass 
splitting pattern of the p and A\ states we are seeing in the chiral limit as Nf is varied. 

Our findings at Nf = 12 are in disagreement with [32, 33]. Lessons from the Dirac spectra and 
RMT to complement p-regime tests are discussed in the next section. 

7. Epsilon regime, Dirac spectrum and RMT 

If the bare parameters of a gauge theory are tuned to the £ -regime in the chirally broken phase, 
the low-lying Dirac spectrum follows the predictions of random matrix theory. The corresponding 
random matrix model is only sensitive to the pattern of chiral symmetry breaking, the topological 
charge and the rescaled fermion mass once the eigenvalues are also rescaled by the same factor 
~L con dV . This idea has been confirmed in various settings both in quenched and fully dynamical 
simulations. The same method is applied here to nearly conformal gauge models. 

The connection between the eigenvalues A of the Dirac operator and chiral symmetry breaking 
is given in the Banks-Casher relation [71], 

^cond = = lim Mm Mm 

A->0m— >0V— ><=o V 
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Figure 7: From simulations at Nf = 4 the first column shows the approach to quartet degeneracy of the 
spectrum as j3 increases. The second column shows the integrated distribution of the two lowest quartets 
averaged. The solid line compares this procedure to RMT with Nf = 4. 

where ~L C ond designates the quark condensate normalized to a single flavor. To generate a non- 
zero density p(0), the smallest eigenvalues must become densely packed as the volume increases, 
with an eigenvalue spacing AA wl/p(0) = 7t /(L conc iV). This allows a crude estimate of the quark 
condensate L com /. One can do better by exploring the £ -regime: If chiral symmetry is spontaneously 
broken, tune the volume and quark mass such that 4- -C L -C A- , so that the Goldstone pion is 
much lighter than the physical value, and finite volume effects are dominant as we discussed in 
Section 2. The chiral Lagrangian of Eq. (2.2) is dominated by the zero-momentum mode from 
the mass term and all kinetic terms are suppressed. In this limit, the distributions of the lowest 
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eigenvalues are identical to those of random matrix theory, a theory of large matrices obeying 
certain symmetries [72, 73, 74]. To connect with RMT, the eigenvalues and quark mass are rescaled 
as z = AZ com /V and pL = m q Jl com iV, and the eigenvalue distributions also depend on the topological 
charge v and the number of quark flavors Nf. RMT is a very useful tool to calculate analytically all 
of the eigenvalue distributions [75]. The eigenvalue distributions in various topological sectors are 
measured via lattice simulations, and via comparison with RMT, the value of the condensate L com i 
can be extracted. 

After we generate large thermalized ensembles, we calculate the lowest twenty eigenvalues 
of the Dirac operator using the PRIMME package [76]. In the continuum limit, the staggered 
eigenvalues form degenerate quartets, with restored taste symmetry. The first column of Fig. 7 
shows the change in the eigenvalue structure for Nf = 4 as the coupling constant is varied. At j8 = 
3.6 grouping into quartets is not seen, the Goldstone pions are somewhat still split, and staggered 
perturbation theory is just beginning to kick in. At j3 = 3.8 doublet pairing appears and at j8 = 4.0 
the quartets are nearly degenerate. The Dirac spectrum is collapsed as required by the Banks- 
Casher relation. In the second column we show the integrated distributions of the two lowest 
eigenvalue quartet averages, 



JO 

which is only justified close to quartet degeneracy. All low eigenvalues are selected with zero 
topology. To compare with RMT, we vary = m ? E con ^y until we satisfy 



where (Ai) smi is the lowest quartet average from simulations and the RMT average (z)rmt depends 
implicitly on jU and Nf. With this optimal value of jU, we can predict the shapes of Pk{X) and 
their integrated distributions, and compare to the simulations. The agreement with the two lowest 
integrated RMT eigenvalue shapes is excellent for the larger /3 values. 

The main qualitative features of the RMT spectrum are very similar in our Nf = 8 simulations 
as shown in Fig. 8. One marked quantitative difference is a noticeable slowdown in response to 
change in the coupling constant. As j8 grows the recovery of the quartet degeneracy is considerably 
delayed in comparison with the onset of p-regime Goldstone dynamics. Overall, for the Nf = 4,8 
models we find consistency between the p-regime analysis and the RMT tests. Earlier, using Asqtad 
fermions at a particular j3 value, we found agreement with RMT even at Nf = 12 which indicated a 
chirally broken phase [24]. Strong taste breaking with Asqtad fermions leaves the quartet averaging 
in question and the bulk pronounced crossover of the Asqtad action as /3 grows is also an issue. 
Currently we are investigating the RMT picture for Nf = 9, 10, 11, 12 with our much improved 
action with stout smearing. This action shows no artifact transitions and handles taste breaking 
much more effectively. Firm conclusions on the Nf = 12 model to support our findings of ^SB in 
the p-regime will require continued investigations. 

8. Inside the conformal window 

We start our investigation and simulations of the conformal window at Nf = 16 which is the 
most accessible for analytic methods. We are particularly interested in the qualitative behavior of 
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Figure 8: The solid lines compare the integrated distribution of the two lowest quartet averages to RMT 
predictions with Nf = 8. 

the finite volume spectrum of the model and the running coupling with its associated beta function 
which is expected to have a weak coupling fixed point around g* 2 0.5, as estimated from the 
scheme-independent, two-loop beta function [77]. A distinguished feature of the Nf = 16 confor- 
mal model is how the renormalized coupling g 2 {L) runs with L, the linear size of the spatial volume 
in a Hamiltonian or Transfer Matrix description. On very small scales the running coupling g 2 {L) 
grows with L as in any other asymptotically free theory. However, g 2 (L) will not grow large, and 
in the L — > oo limit it will converge to the fixed point g* 2 which is rather weak, within the reach 
of perturbation theory. There is non-trivial, small- volume dynamics which is illustrated first in the 
pure gauge sector. 

At small g 2 , without fermions, the zero-momentum components of the gauge field are known 
to dominate the dynamics [78, 79, 80]. With SU(3) gauge group, there are twenty-seven degen- 
erate vacuum states, separated by energy barriers which are generated by the integrated effects 
of the non-zero momentum components of the gauge field in the Born-Oppenheimer approxima- 
tion. The lowest-energy excitations of the gauge field Hamiltonian scale as ~ g 2 / 3 (L) jh evolving 
into glueball states and becoming independent of the volume as the coupling constant grows with 
L. Non-trivial dynamics evolves through three stages as L grows. In the first regime, in very 
small boxes, tunneling is suppressed between vacua which remain isolated. In the second regime, 
for larger L, tunneling sets in and electric flux states will not be exponentially suppressed. Both 
regimes represent small worlds with zero-momentum spectra separated from higher momentum 
modes of the theory with energies on the scale of 2n/L. At large enough L the gauge dynamics 
overcomes the energy barrier, and wave functions spread over the vacuum valley. This third regime 
is the crossover to confinement where the electric fluxes collapse into thin string states wrapping 
around the box. 

It is likely that a conformal theory with a weak coupling fixed point at Nf = 16 will have only 
the first two regimes which are common with QCD. Now the calculations have to include fermion 
loops [81, 82]. The vacuum structure in small enough volumes, for which the wave functional 
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is sufficiently localized around the vacuum configuration, remains calculable by adding in one- 
loop order the quantum effects of the fermion field fluctuations. The spatially constant abelian 
gauge fields parametrizing the vacuum valley are given by A,-(x) = T a Cf /L where T a are the (N-l) 
generators for the Cartan subalgebra of SU{N). For 5*7(3), T\ = X 3 /2 and T 2 = Ag/2. With Nf 
flavors of massless fermion fields the effective potential of the constant mode is given by 

V e k ff(C fe ) = £V(C fe [ M « -liPV-NfYyl&tf + xk), (8.1) 

i>j i 

with k = for periodic, or k = (1, 1, 1), for antiperiodic boundary conditions on the fermion fields. 
The function V(C) is the one-loop effective potential for Nf = and the weight vectors juW are 
determined by the eigenvalues of the abelian generators. For SU(3) jJ.^ = (1, 1 , — 2) / \/T2 and 
= 2 ( 1 , — 1 , 0) . The correct quantum vacuum is found at the minimum of this effective potential 
which is dramatically changed by the fermion loop contributions. The Polyakov loop observables 




Figure 9: The time evolution of complex Polyakov loop distributions are shown from our Nf = 16 simula- 
tions with 12 3 x 36 lattice volume. Tree-level Symanzik-improved gauge action is used in the simulations 
and staggered fermions with three stout steps and very small fermion masses. 

remain center elements at the new vacuum configurations with complex values; for SU (N) 

Pj = Itr (exp(/C^)) = ^ J>p(irf°c5) = exp(2Kilj/N). (8.2) 

This implies $ } C b = 2%\/N (mod 2k), independent of n, and V e ^ = -N f NV(2n\/N + rik). In the 
case of antiperiodic boundary conditions, k = (1,1,1), this is minimal only when 1 = (mod 2%). 
The quantum vacuum in this case is the naive one, A = (Pj = 1). In the case of periodic boundary 
conditions, k = 0, the vacua have 1 ^ 0, so that Pj correspond to non-trivial center elements. For 
SU(3), there are now 8 degenerate vacua characterized by eight different Polyakov loops, Pj = 
exp(±27T//3). Since they are related by coordinate reflections, in a small volume parity (P) and 
charge conjugation (C) are spontaneously broken, although CP is still a good symmetry [81]. 

Our simulations of the Nf = 16 model below the conformal fixed point g* 2 confirm the theo- 
retical vacuum structure. Fig. 9 shows the time evolution of Polyakov loop distributions monitored 
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along the three separate spatial directions. On the left side, with periodic spatial boundary condi- 
tions, the time evolution is shown starting from randomized gauge configuration with the Polyakov 
loop at the origin. The system evolves into one of the eight degenerate vacua selected by the 
positive imaginary part of the complex Polyakov loop along the x and y direction and negative 
imaginary part along the z direction. On the right, with antiperiodic spatial boundary conditions, 
the vacuum is unique and trivial with real Polyakov loop in all three directions. The time evolution 
is particularly interesting in the z direction with a swing first from the randomized gauge configu- 
ration to a complex metastable minimum first, and eventually tunneling back to the trivial vacuum 
and staying there, as expected. The measured fermion-antifermion spectra and the spectrum of the 
Dirac operator further confirm this vacuum structure. Our plans include the continued investiga- 
tion of zero-mode gauge dynamics which should clarify many important aspects of conformal and 
nearly conformal gauge theories. 
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